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ABSTRACT 



We analyze the window functions for the spherical harmonic mode estimators of all-sky, 
volume limited surveys considering evolutionary effects along the past light-cone which include 
the deviation of the distance scale from a linear relationship with rcdshift, linear peculiar 
velocity corrections, and linear evolution of the density perturbations. The spherical harmonic 
basis functions are considered because they correspond most closely to the symmetries of typical 
survey geometries and of the light-cone effects we consider. Our results show substantial 
broadening of the windows over that expected by ignoring light-cone effects, indicating the 
difficulty of measuring the power spectrum independently from cosmology. We suggest that 
because of light-cone effects, deep redshift surveys should either be analyzed in conjunction 
with CMBR data which determines the cosmological parameters, or by using a Bayesian 
likelihood scheme in which varying cosmological parameters and a simple parameterization of 
the primordial power spectrum are assumed as the priors, so that observed data can be mapped 
from redshift to real space. The derived power spectrum can then be compared to underlying 
models of fluctuation generation and growth in structure formation to evaluate both these 
models and the cosmological priors. 
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The ability of a number of current and proposed galaxy surveys, e.g. The Sloan Digital Sky Survey 
(hereafter SDSS, Gunn fc Weinberg (1995), Strauss (1997)[ ) and the two-degree-field survey (Taylor & Gray 



(1994)), to probe deep redshifts offers exciting opportunities to observe structure on the largest scales, 



yielding results which will be invaluable for cosmology. As is often the case, however, the data collected 
is not directly compatible with that which theorists calculate. To see this, consider a typical theoretical 
calculation of large scale structure. Structure is presumed to have evolved from the gravitational collapse of 
small initial fluctuations, and on scales larger than approximately 10 Mpc, this collapse is well described by 
linear perturbation theory. Usually, theorists express the outcomes of their structure-formation models in 
terms of the power spectrum, the Fourier transforms of the two-point correlation function given in conformal 
spatial coordinates, evaluated on a constant time slice. That is, they measure the power spectrum of the 
structure for a snapshot of the universe taken at a fixed time. A galaxy survey, however, detects structure 
on a light-cone time slice, and measures redshifts, not conformal distances. For lengths much shorter than 
the Hubble length, the distinction is unimportant since to first order distance is proportional to redshift, 
and structure does not evolve much over the light-crossing time of the survey. But, the newest surveys, 
such as the SDSS, will penetrate to as much as 50% of the Hubble distance. Either assuming that proper 
distance is linearly proportional to redshift over these distance scales, or ignoring the evolution of structure 
over the corresponding time scales, leads to significant errors in the estimation of the power spectrum. In 
this letter, we examine these effects and determine their significance for estimating the power spectrum 
from a deep galaxy survey. 

All surveys are limited in the volume of space which they can cover, and so are less than ideal 
observations. Furthermore, shot noise arises from the finite sampling inherent in any galaxy survey. 
Conversely, theoretical predictions for power spectra are determined for the whole universe continuously 
sampled and are therefore expressed in terms of a continuous power spectrum (unless the universe is 
spatially compact, but no current evidence suggests that it is). In sampling a finite volume, one probes 
instead a discrete set of eigenmodes and the amplitudes of the observed modes will be a convolution of the 
continuum modes with a window function. Ideally, this window function is narrow so that the observed 
mode amplitudes correspond closely to the their continuum counterparts, and this can be achieved through 
the proper weighting of the observational data. 



For Fourier decomposition, Fegmark (1995) has developed a method of analyzing redshift surveys 
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guaranteed for a specified survey geometry to minimize the width of the window function to the shot noise 
limit. Fourier decomposition, however, is not well suited to the analysis of most redshift surveys because 
survey geometries are usually defined by slices of a sphere, and because redshift effects depend only on the 
distance from the observer. A better basis is one which takes advantage of this symmetry, i.e. the spherical 



harmonic basis ( Heavens and Talor (1995) ). For standard homogeneous, isotropic models, the power 
spectrum of the spherical harmonic decomposition is equivalent to that of the Fourier decomposition — 
each depends only on the magnitude of the mode considered, so theoretical calculations are equally useful 
for each basis. Our analysis therefore addresses the determination or estimation of the power spectrum in 
the spherical harmonic basis. 

Given a survey geometry, we would like to find a weighting function ip which will produce the "optimal" 
estimate for the amplitude of a given £q, too, and ko harmonic of the continuum fluctuation spectrum. By 
optimal we mean that ip should simultaneously minimize the width of the window function in £, to, and k 
space and minimize the shot noise signal due to finite sampling — under the constraint that ip is non-trivial. 



Following Tegmark (1995), we define the penalty function 

IP - PI 2 

|l-4| 2 + |m-mol 2 + |fc 2 fc °' , (1) 

7 

which forces our window toward the £o,mo, and fco we are trying to observe. The shot noise is proportional 
to n(z) -1 , where n(z) is the average number of galaxies observed in a unit volume at redshift z. To find the 
optimal ip we solve the eigen value problem 

(\£ - £ \ 2 + \m- too| 2 + |fc2 ' 2 fc ° 2 ' 2 + V = (2) 

where the ground state represents the optimal weighting function as it minimizes E. We normalize it, for 
convenience, by requiring 

' d 3 z^ 2 (z) = 1. (3) 



Note that tjj can always be chosen to be a real valued function (which we shall assume it is), simplifying 
the numerics. Also note that for the radial portion of the penalty function we choose |P — A;q | 2 instead 
of \k — fco| 2 , which one might naively select. We do so because there is no simple way to express k in 
coordinate space; the operator P is simply negative the Laplacian in physical space. The operators £ and to 
return the magnitude of the total angular momentum and the z component respectively. Since the angular 
portion is separable, we can re-express ip as a product ipiiuk(z)Y tJ-l/ (Q z ) where /i and v are the eigenvalues of 
the generalized spherical harmonics which satisfy the boundary conditions given by the survey geometry. 
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In the limit of an all sky survey, ji = £o arL d v — tuq which are integers, although this not generally so. 
The parameter 7 controls how we wish to weight shot noise to the window width. By increasing 7, we can 
reduce shot noise by accepting a larger variation in k , or by reducing 7 we can get a narrower k variation, 
but at the expense of increased shot noise. For the special case of a volume limited survey in which n(z) 
is a constant, the shot noise can be absorbed into the right hand side of the equation and the choice of 7 
is irrelevant; there is only one optimal ip. As we have already pointed out, the angular part of the solution 
is given by the generalized spherical harmonics. The solution to the radial portion is then just a spherical 
Bessel function j^(k^z), where j^{k^z max ) — at the outer boundary of the volume z max , and k^ is the zero 
closest to k, the wavenumber of the mode whose amplitude we wish to estimate. It is this simplified problem 
which we will consider in this letter since the numerics are easier, and we expect that any light-cone effects 
which we observe in a more general survey of similar depth will be at least if not more severe. 

Given the optimal weighting function, our estimator for the fluctuation amplitude of the k, I, m mode 

is 

a lm (k) = [ ^(^-n^ 
J n(z) 

where n{z) = (n(z)} is the ensemble average of the number density of galaxies at redshift z. A more 
interesting quantity is the square of this estimator, which averaged over the ensemble of realizations should 



be related to the continuum power spectrum. Using the analysis of [Peebles (1980)| , one can straight 
forwardly show that in the ensemble average 

(&lm(k)km(k)) = [ d 3 z'^- + [ d 3 zd 3 z'^z,z')^(z)^(z'), (5) 
J n{z) J 

where £ is the correlation function in redshift space. From the above, we can infer an estimator for the 
power in a particular mode 

P im (k) = a lm (k)d lm (k) - f d 3 z^$, (6) 



n(z) 

where we subtract out the expected shot noise contribution. The light cone projection of the correlation 
function for flat space (hyperbolic, or "open" models will be covered in a more detailed paper |de Laix"fc 



Starkman (1998) ) has been calculated by Nakamura, Matsubara, fc Suto (1997) . In terms of a spherical 
harmonic basis, £ may be expressed as 

/2k 2 dk 
—^D(z)D(z')Y em (n z )Y em (n' z ) (7) 

: . . 



Im 



H(z) ( k£ . ... fc(l+l) . 
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k 2 {A e jt-2(kr) - Bij e (kr) + C e jt- 2 {kr)) 



with 



B t 
C e 



£{£-!) 



(2^ + 1)(2£ - 1) 

e 

(2€ + 1)(2£ - 1) 
(I + !)(£ + 2) 
(2£ + l)(2£ + 3)' 



(2i+l)(2£+3) 



The proper radial distance r is 



dz 



io H(z) 

where we introduce the redshift-dependent Hubble parameter for flat space 



H{z) = H y/n (l + zf + A . 



(8) 



(9) 



The parameters f2o and Ao represent respectively the present day matter density and present day 
cosmological constant, in units of the critical density, and f^o + Ao = 1 for flat space. We also use the linear 
growth factor 



= 5 -^l H (z) r ir^dz', (io) 



H(z') 



along with the linearly evolving bias parameter which, following Fry (1996) , can be written 



h{z)=i+ W) {b °~ 11 



(11) 



to define an overall linear growth factor for the galaxy perturbations relative to the underlying density 
perturbations: 



D(z) 



D(z)b(z) 



D(o)b(oy 

Finally, we use linear velocity distortion parameter (3 given by 

d\nD(z) 1 



(12) 



(3{z) 



d\n{l + z) b(z)' 

The reader is reminded that eq. (0) is only valid when Qo + Ao = 1, otherwise an expression with the 
generalized spherical Bcssel functions derived from a hyperbolic Laplacian are required. 

Using our expression for the correlation function, we can show that the ensemble average of the 
estimator in eq. (^|) is equal to 



(13) 



p lm {k))= { fc' 2 dfc'p(fc')|E 7 '™'( fc ')' 



(14) 
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with 



I lm {k) = / d 3 z^{z)D(z)Y em {n z ) (15) 



J • a. ^ 



H{z) ( kl . fc(l+l) . 

U-i(kr) - — — — 31+iikr) 



1 + z \2£+V K ' 21+1 

- k 2 (A e je-2(kr) - B t j t (kr) + Gt7/_ 2 (Ar)) 

In the special case of an all sky survey, where ip(z) oc Yi omo (Q z ) 1 one can perform the angular integral dfl z 
trivially. The estimator for the power in a particular mode is then reduced to the convolution of the true 
continuum power spectrum with a window function 

W ek = fc 2 /|(fc), (16) 

which is independent of m. It is these windows which indicate how well a power spectrum derived from a 
limited rcdshift survey will correspond to the underlying continuum power spectrum. 

To demonstrate how light-cone effects will influence the observational window functions defined above, 
we compare the windows which arise when all redshift effects are included to those obtained ignoring 
light-cone effects by taking the distance r proportional to the redshift (r = H Q ' 1 z) and holding the other 
redshift dependent parameters — /3, H, and D - fixed at their z = values. Often one sees results calculated 
in the latter naive limit. For simplicity we consider only volume-limited all-sky surveys, investigating two 
survey depths: z max = 0.25 and z max = 0.5, which correspond to the approximate depths of the SDSS 
magnitude- limited galaxy survey and the SDSS bright red galaxy survey. We take unbiased (bo = 1), Qq = 1 
cosmologies, and consider two values of k, 0.033 and .33 h Mpc -1 . 

In figure ^, we show a plot of the naive window functions for k — 0.033 h Mpc" 1 with a range of £'s in 
a z < 0.25 volume-limited survey, ignoring the light-cone effects, i.e. this shows the optimal results one can 
hope to achieve for a volume-limited survey of the specified depth. In figure ^|, we show the same results for 
a realistic survey, including light -cone effects. In figures | and |, we consider the same survey geometry, but 
calculate the window functions for k — 0.33 h Mpc^ 1 . The former plot shown the naive window function 
while the latter shows the realistic result. Next we examine the deeper survey geometry which extends to 
a depth of z = 0.5. Again we consider first the larger scale, k = 0.033 h Mpc" 1 , in the naive and realistic 
limit in figures |^ and || respectively. Finally, we examine the smaller scale in figures ^ and |^. 

The results shown in the figures described above reveal the importance of considering light-cone effects 
when estimating power from galaxy surveys, with several effects meriting closer attention. The most striking 
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is the broadening of the realistic window functions compared to their idealized counterparts. This is caused 
predominantly by the deviation from the linear distance-redshift relation. At higher rcdshifts, the larger 
Hubble expansion rate causes the galaxies to appear more spread out than they do locally. Thus a wave that 
probes a given scale locally, probes a smaller scale (larger k) at high redshift. A mode in conformal space 
will then couple to a range of scales when mapped into redshift space, which leads to a broadening in the 
window. It also causes a readily observed shift in the peak to larger scales. We can estimate the broadening 
quite easily if we consider the effective frequency at which a conformal mode oscillates in redshift space. For 
large values of the argument, a spherical Bcssel function can be approximated ji{kr) ks sin(fcr + £tt /2)/kr. 
In a small region Az located at z, we can rewrite this as sin(/cr(z) + kdr/dzAz + £ir/2)/kr(z), which has an 
effective frequency of oscillation of kdr/dz. Recall from eq. (JsJ) that dr/dz = H^ 1 (z). This tells us that the 
effective frequency of the conformal modes which describe the density fluctuations will drift from an initial 
value of k to H~ l (z max )k in redshift space, so the inverse of this quantity should be a good estimate of the 
actual window width. For our examples of z max = 0.25 and 0.5, we see that the window widths should be 
about 1.4fc and 1.8k respectively, which are in good agreement with the results presented in the figures. A 
final effect, which is less visible, is an overall decrease in the amplitude of the window function. Since the 
power spectrum decreases with growing z, the overall power measured in a given mode will be less than its 
local value. 

These results are important for analyzing observations, because they indicate that it is impossible to 
make a clean measurement of the power spectrum, independent of cosmology, from a deep redshift survey. 
The light-cone effects result in significant broadening and shifting of the window functions when compared 
to their idealized limit, in a way that depends sensitively on the cosmology. One must deconvolve these 
effects to compare observations to theory, but this requires knowledge of the underlying cosmological model 
(perhaps soon to be provided by cosmic microwave background observations). An alternative approach 
might be to do a Bayesian analysis assuming different cosmologies as priors. One maps the observed data 
into conformal coordinate space, utilizing the selected cosmology, calculates a power spectrum estimate, 
and compares the estimate to theoretical predictions based on simple parameterizations of the underlying 
power spectrum from theories like inflation or defects. For the correct cosmological model, the windows 
will have minimum width, while the incorrect models will have broadened windows. It is possible that 
these effects can be used to learn about the underlying cosmology, and even to break the degeneracy in the 
power spectrum for models in which the shape parameter, Vlh, is constant (we are currently investigating 
this problem). 
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Figure Captions 

Fig. 1. — The naive window function (no light-cone effects) for an all-sky survey of depth z = .25 used to 
optimally estimate the power in the .033 h Mpc -1 mode. The results for number of different I modes are 
shown. 

Fig. 2. — The true, light-cone corrected window functions for an all-sky survey of depth z = .25 used to 
optimally estimate the power in the .033 h Mpc -1 mode. The results for number of different I modes are 
shown. 

Fig. 3. — The naive window function (no light-cone effects) for an all-sky survey of depth z = .25 used to 
optimally estimate the power in the .33 h Mpc -1 mode. The results for number of different £ modes arc 
shown. 

Fig. 4. — The true, light-cone corrected window functions for an all-sky survey of depth z = .25 used to 
optimally estimate the power in the .33 h Mpc -1 mode. The results for number of different I modes arc 
shown. 

Fig. 5. — The naive window function (no light-cone effects) for an all-sky survey of depth z = .5 used to 
optimally estimate the power in the .033 h Mpc -1 mode. The results for number of different I modes are 
shown. 

Fig. 6. — The true, light-cone corrected window functions for an all-sky survey of depth z = .5 used to 
optimally estimate the power in the .033 h Mpc -1 mode. The results for number of different £ modes are 
shown. 

Fig. 7. — The naive window function (no light-cone effects) for an all-sky survey of depth z = .5 used to 
optimally estimate the power in the .33 h Mpc -1 mode. The results for number of different £ modes are 
shown. 

Fig. 8. — The true, light-cone corrected window functions for an all-sky survey of depth z = .5 used to 
optimally estimate the power in the .33 h Mpc -1 mode. The results for number of different £ modes are 
shown. 



